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Abstract 

We study the asymptotic zero distribution of type II multiple orthogonal poly- 
nomials associated with two Macdonald functions (modified Bessel functions of the 
second kind). Based on the four-term recurrence relation, it is shown that, after 
proper scaling, the sequence of normalized zero counting measures converges weakly 
to the first component of a vector of two measures which satisfies a vector equilib- 
rium problem with two external fields. We also give the explicit formula for the 
equilibrium vector in terms of solutions of an algebraic equation. 

1 Introduction 

Given a positive measure on the real line for which the support is not finite and all 
moments exist, there exists a sequence of monic orthogonal polynomials {pk} of degree 
k. Such polynomials satisfy a three-term recurrence relation of the form 

xpk{x) = pk+i{x) + bkPk{x) + alpk-i{x), k>0, (1.1) 

with afc > 0, 6a; G M and po = 1, p-i = 0. 

It is well-known that the zeros of are real and simple. We can associate with Pk{x) 
the normalized zero counting measure 



1 ^ 

^Pk) = (1.2) 



where Xj^k, j = 1, . . . ,k, are the zeros of pk and 5x denotes the Dirac point mass at x. 
A measure v is called the asymptotic zero distribution of {pk} if 



lim / f du{pk) = I fdv 
k^coj J 
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for every bounded continuous function / on M, i.e., it is the weak limit of the measures 

Suppose hmfc_^oo o,k = a and Hm^^oo bk = b with a > and 6 G M, then the polyno- 
mials generated by (jl.ip have the asymptotic zero distribution ?^[a,/3] with density 

^^["'^]' (1.3) 

^•'^ 1 0, elsewhere, 

where a = b — 2a, (3 = b + 2a; cf. [16] . 

This result has been extended in |15) to the case of orthogonal polynomials generated 
by the recurrence relation 

XPk,n{^) = Pk+l,ni^) + bk,nPk,n{x) + al „Pk-l,n{x) , k,nGn, 

with varying recurrence coefficients ak^n > and bk.n G ^ depending on a parameter n. 
Assume that the recurrence coefficients have continuous limits 

lim ak,n = a{s), lim bk,n = b{s), 

k/n^s k/n^s 

where a : (0,oo) — t- [0,oo), b : (0,oo) — t- M and the notation liuik/n^s means that both 
A;,n — )■ oo with k/n ^ s > 0, it is proved that the asymptotic zero distribution is given 
by the average 

1 

lim u{pk-n) = - / W[a{s),(Sis)]ds, (1.4) 
k/n—^s S Jq 

where a{s) := b{s) — 2a{s), /3{s) := b{s) + 2o(s) and 'W[a,i3] is defined by (|1.3p if a < /3 
and hy 5a if a = /3; see Theorem 1.10 of [15]. 

A natural generalization of this case consists in considering, for each n G N, polyno- 
mials Pk^n satisfying an m-term recurrence relation with varying coefficients: 

xPk,nix) = Pk+lAx) + btiPk,n{x) + b^kiPk-lAx) + ■■■+ 6i";"'^Pfc+2-m,n(x), (1.5) 



where Pq = 1, P_i = 0, ...,P_m+2 = and the recurrence coefficients have scaling 
limits 

fc/n— >s 



lim bf^ = b^^\s), j = 0,...,m 



for certain functions 6*-'^^ , • • • , fe^™ • 

For the simplest case b^^^ = b^^\s), i.e., we remove the dependence on the parameter 
n and the recurrence coefficients in (jl.Sp are actually constant, the zeros of Pk = Pk,n 
are closely related to the spectrum of certain banded Toeplitz matrix. Indeed, if we 
associate with the functions b'^^'^ a family of functions 

As{z) = z + b'^^\s) + 6(^)(s)z-i + • • • + 6("^-2)(s)z-'^+2, (1.6) 
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and the sequence oi k x k Toeplitz matrices {Tj.{As))k with symbol Ag, defined by 



{Tk{As))ji 



1, if/ = j + l, 

6W(s), ifZ = j-«, z = 0,...,m-2, (1.7) 
0, otherwise, 



it is readily seen that Pk{^) = if and only if A is an eigenvalue of Tk{As). Hence, the 
investigation of limiting zero distribution of is equivalent to the study of the limiting 
behavior of the spectrum of Tk{As) as — >■ oo. 

The limiting behavior of the spectrum of Tj^{As) as A; — )• oo is characterized by the 
solutions of the algebraic equation As{z) = x; see [3]. For every x G C, there exist 
exactly m — 1 solutions of the equation As{z) = x (assume that ?)(™~^)(s) 7^ 0), which 
we denote hy Zj{x, s), j = 1, ... ,m — 1 and label these solutions by their absolute value 
so that 

\zi{x,s)\ > 1^2(3;, s)| > • • • > \zm-iix,s)\ > 0. (1.8) 

We put 

Ti{s) = {xeC \ \zi{x, s)\ = \z2{x, s)\}, (1.9) 

which is a finite union of analytic arcs. 

It was shown by Schmidt and Spitzer [20] that the eigenvalues of Tk{As) accumulate 
on the contour Ti{s) as k tends to 00. Moreover, Hirschman [12J proved that the sequence 
of normalized counting measures of the eigenvalues of (Ag) converges weakly to a Borel 
probability measure /xf supported on Ti{,s) as A; — )• 00; see also |31 Chapter 11]. The 
precise form of /if is given by 

d,Ux) = ^ ('-^ - dx, (1.10) 

'^^^ ^ 2Tri \zi-{x,s) zi+{x,s)J ^ ' 



which is due to the result in [TO]. Here, we have that ' denotes the derivative with respect 
to x, dx is the complex line element on ri(s) and z\±{x^ s) is the limiting value of z\{x^ s) 
as X — )• X from the it side of ri(s). Moreover, the measure /if is also characterized by 
an equilibrium problem; see Theorem 12.21 below for a statement in the context of the 
specific example considered in this paper. 

Under certain conditions, the polynomials Pfc,n satisfying the recurrence (jl.Sp have 
a limiting zero distribution as well, which is an average, with respect to the parameter 
s, of the measures (|1.1U|) . More precisely, we have (see Theorem 1.2 in [,14j ) : 

Theorem 1.1. Let for each n G N, m — 1 sequences {b^j^l^}'^^, j = 0, . . . , m — 2, of real 

coefficients be given and assume that there exist continuous functions b^^^ : [0, 00) — t- R, 
j = 0, . . . ,m — 2, such that for each s > 0, 

lim bl^l = b^^\s), j = 0,...,m-2. (1.11) 

k/n^s ' 

Let Pk n be the monic polynomials generated by the recurrence (jl.Sp and suppose that 
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(a) 



The polynomials Pk,n have real and simple zeros x^"" < ■ ■ ■ < x^"" satisfying for 
each k and n the interlacing property 



k+l,n 
'j 



. k,n , 
< Xj < X 



k+l,n 



for j = 1, . . . 



k 



(b) Ti(s) C M for every s > 0, where Ti{s) is given by ()1.9p . 

Then the normalized zero counting measures v{Pk,n) = \ Z]y=i ^ have a weak limit 
as k,n ^ oo with /c/n — )■ ^ > given by 



where fif is the measure (jl.lOp . 

It is worth noting that in fl], the authors present a conditional theorem giving the 
asymptotic zero distribution for polynomials satisfying a specific four-term recurrence 
relation. 

The aim of this paper is to give more insight on the nature of the asymptotic zero 
distribution in the particular case of type II multiple orthogonal polynomials associated 
with two Macdonald functions (modified Bessel functions of the second kind) K^{x) 
{v > 0). A feature of the present case is the appearance of a vector equilibrium problem 
with two external fields, for which the first component of the unique minimizer is the 
weak limit of the normalized counting zero measures. 

We mainly follow the idea in |14j . where the authors consider a model of non- 
intersecting squared Bessel paths and derive a vector equilibrium problem for the limiting 
zero distribution of type II multiple orthogonal polynomials associated with the modi- 
fied Bessel functions of the first kind [6ll3- In that case, the vector equilibrium problem 
involves two measures supported on the positive real line and the negative real line, 
respectively, with an external field acting on the first measure and a constraint acting 
on the second measure |141 Theorem 1.7]. 

2 Statement of results 

2.1 Multiple orthogonal polynomials associated with Macdonald func- 




(1.12) 



tions 



Assuming x > 0, we define the scaled Macdonald function pi, by 



p^{x) = 2x''/^K^{2,/^) 



(2.1) 



and consider two weights 



dfii{x) = x°' pi,{x)dx, dfi2{x) = x'^pu+i{x)dx, a > —1, > 0, 



(2.2) 
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on the positive real line. For any fc, m G N, the type II multiple orthogonal polynomials 
Pkm system of weights (/^i, /X2) are such that is a monic polynomial of degree 

k + ra and satisfies the following multiple orthogonality conditions: 



= 0, i = 0,l,...,A:-l, (2.3) 

) 

p^^„(x)x-'d/U2(x) = 0, J =0,l,...,m- 1. (2.4) 

By taking m = k, we set 

P2kix) = Pt,k{x), P2k+i{x) = Pt+i,kix). 
An explicit formula for is given by 

k 

Pk{x)=^ak{j)x'-^, (2.5) 
j=0 

where 

'k\ {a + l)kia + u + l)k 



\jj {a + l)k-j{a + + l)k-j 

see O Theorem 2] . It is shown in [23] that Pk satisfies the following four-term recurrence 
relation 

xPkix) = Pk+i{x) + bkPk{x) + CkPk-iix) + dkPk-2{x) (2.6) 
with recurrence coefficients 

hk = {k + a + l){?>k + a + 2v) -{a + l){v- 1), 

c,, = k{k + a){k + a + v){?,k + 2a + u), (2.7) 
dj^ = k{k - l){k + a - l){k + a){k + a + V - l){k + a + u). 

These polynomials constitute one of few examples of multiple orthogonal polynomials 
that are not related to the classical orthogonal polynomials. They are first introduced 
by Van Assche and Yakubovich in [23j, which solve an open problem posed by Prunikov 
[18] ; see also [21 H] for recent study. 

Our goal is to investigate the limiting zero distribution of scaled polynomials Pk- 
Namely, we introduce a new parameter n E N and put 

P^nix) := (2.8) 

Clearly, Pk,n{x) is a polynomial of degree k for each n. In view of (j2.6p ~ (j2.8|) . it is readily 
seen that Pk,n{x) satisfies the following recurrence relation 

xPk,n{x) = Pk+l,n{x) + bk,nPk,n{x) + Ck,nPk-l,n{^) + dk,nPk-2,ni^) , (2.9) 
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with recurrence coefficients given by 

_ {k + a + l){3k + a + 2u)-{a + - 1) 

k(k + a)(k + a + u)(3k + 2a + i') /r,1^^ 

Ck,n = 4 , {^■^^) 

_ k{k-l){k + a-l){k + a){k + a + v -l){k + a + v) 
As in p.2p . the normalized counting zero measure of Pk,n is defined by 

^{Pk,n) = \ Yl (2.11) 

We will derive a vector equilibrium problem with two external fields and show that the 
first component of the equilibrium vector is the weak limit of y{Pk,n) as A;, n — >• oo with 
fc/n — )• > 0. The equilibrium vector itself can be explicitly given in terms of the 
solutions of an algebraic equation. Our results are actually rather general in the sense 
that we also allow the parameters a and u to increase proportionally to n as n increases. 

2.2 Statement of results 

We scale the parameters a and v in the following way: 

a pn, V ^ qn, (2-12) 

with p,q > 0. The results corresponding to a and i' fixed can be obtained by taking the 
limits as p, g — >• 0, respectively. 

Let A;, n — 7- oo in such a way that k/n — >• s, for some s > 0, we then observe that the 
recurrence coefficients ()2.10p have scaling limits b{s), c{s) and d{s) given by 

lim bk,n = b{s) = 3s^ + Asp + 2sq + p'^ + pq, 

k/n^s 

lim Ck,n = c{s) = s{s+p){s+p + q){3s + 2p + q), (21^) 
lim = = s^(s+p)^(s+p + g)^. 

Clearly, these limits depend on p and q. As in ()1.6p . we have the associated family 
of symbols 

As{z) = z + b{s) + c{s)z-^ + d{s)z-\ (2.14) 

and the solutions zi{x,s), Z2{x,s) and 23(3;, s) of the algebraic equation As{z) = x. We 
define ri(s) as in (|1.9p and similarly 

T2{s) = {x G C I 1^2(3;, s)| = |z3(x, s)\]. (2.15) 

The following proposition ensures that the polynomials Pk,n satisfy the hypothesis 
(a) of Theorem 11.11 
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Proposition 2.1. Let p,q > 0. Then the polynomials Pk,n generated by (|2.9|) with 
recurrence coefficients ()2.10p have real and simple zeros in (0, oo) with the interlacing 
property. 

Proof. The proof of the proposition follows from the fact that the measures ((i/ii,d/X2) 
from ()2.2p form an AT system (cf. [221 S])) which implies all the zeros of Pk are 
simple, lie in (0, +c«) [IS [22], and satisfy the interlacing property [1]. □ 

By Proposition 12.31 stated below, it is easily seen that the hypothesis (b) is also 
satisfied. Therefore, we see from Theorem 1 1 . 1 1 that . with the scaling given in ()2.12p . the 
probability measure 



^l{ds, ^>0, (2.16) 



is the weak limit of the normalized zero counting measures. The main result of this 
paper is that v\ can also be obtained as the first component of a vector of measures 
{u\, i/l) that satisfies a vector equilibrium problem with two external fields. 

To define v^i '^^ need to introduce the second measure /^l, which is supported on 
r2(s) (see ()2.15p ) and given by 

dl^lix) = — - ^+ ]dx, X E r2 s . 2.17 

2m\z2-{x,s) Z2+{x,s)J 

It is a positive measure on r2(s) with total mass 1/2. For each ^ > 0, we define z^l ^ 
manner similar to the definition of in (|2.16p . i.e., 

1 /■« 

;u|ds. (2.18) 



Then i^l is ^ measure on (Js<gr2(s) = r2(^) with total mass 1/2. 

An essential point for the rest of the paper is the main result of pU], which asserts 
that the vector of measures (^f,//2) is characterized by a vector equilibrium problem. 
In the present context, this is stated in the following theorem. 

Theorem 2.2. For each s > 0, the vector (/uf , /n|) is the unique minimizer for the energy 
functional 

\og- — ^ — :dni{x)dfii{y) + II log- — - — :dn2{x)dfi2{y) 



y\ JJ \x-y 

1 



log^ rdfllix)dfi2{y) (2.19) 

\x - y\ 



among all vectors {ij,i,fi2) satisfying supp(/Uj) C ^j{s) for j = 1,2, and 

1 



dm = 1, I dfj.2 2" 
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The measures fif and /i| satisfy the following Euler- Lagrange variational conditions: 
2 J log |x - y\dtiliy) - J log \x - y\diil{y) =t, x G ri(s), (2.20) 
for some constant and 

2 j log\x - y\dfil{y) - j log\x - y\dfiliy) = 0, x£T2is). (2.21) 

We shall obtain the equilibrium problem for the vector of measures (z^^jZ^I) by in- 
tegrating the variational conditions ()2.20p and ()2.2ip with respect to the variable s. 
The main difficulty lies in the fact that Ti{s) and T2{s) are varying with s. Hence, it 
is necessary to study how these contours depend on s. The next proposition reveals 
that Ti(s) and ^2(3) are indeed real intervals and actually are increasing as s increases. 
The monotonicity of Ti{s) and r2(s), as we will see later, plays a important role in the 
derivation of the equilibrium problem. 

Proposition 2.3. For each s > 0, we have that Ti{s) C (0,oo) and T2{s) C (— oo,0). 
More precisely, there exist r/(s) < < /3(s) < 7(5) so that 

r,is) = [/3(s),7(5)], T2is) = {-oo,vis)]. (2.22) 

In addition, we have 

(a) 7(5) is positive and strictly increasing for s > 0, with lims_s.o+ lis) = p{p + q) and 
\mis^cyol{s) = 00, 

(b) (3{s) is positive and strictly decreasing for s > with lims_s.o+ /^(s) = p{p + q) and 
\mis-,oo li{s) = 0, 

(c) "qis) is negative and strictly increasing for s > with lims_^.o+ 7?(s) = —q^ /A and 
lim^^oo ??(s) = 0. 

Figure [U gives an illustrative plot of the functions ^(s), 7(5) and r?(s). 

Now we come to the main result of this paper, i.e., the equilibrium problem for the 
vector of measures (y^^vl). The fact that Ti{s) and r2(s) are increasing as s increases, 
induces two external fields Vi{x) and V2{x) acting on v\ and z/|, respectively. 

Theorem 2.4. For every ^ > 0, the vector of measures (z^^,!/!) is the unique minimizer 
for the energy functional 



II 



log- — ^- — 7dui{x)di>i{y) + [[ log- — - — Tdv2{x)du2{y) 
\x-y\ J J \x-y\ 

log—^Mx)du2{v) + \ I Vi{x)dvi{x) + \ I V2{x)dv2{x), (2.23) 
\x -y\ i J i ' 
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Figure 1: Graph of the curves /3(s), 7(5) and t]{s) with p = 3 and q = 10. 



over all vectors of measures {i'i,i'2) such that supp{i'i) C [0, 00), f dvi = 1 and supp(i/2) C 
(— c«,0], / di'2 = 1/2, where 



Vi{x) = I log 

'0 



2:i(x,s) 



ds and V2 (x) 



/•OD 

/ log 


Z2{x,s) 


ds 


Jo 


Z3{x,s) 





(2.24) 



Z2{X, S) 

The measures and v\ are characterized by the following variational conditions: 



2 J \og\x - y\dvl{y) - j \og\x - y\dvl{y) - ^Vi{x) 



= £, for x G supp(i^^), 
< for X G [0,00), 



(2.25) 



for some i, and 



2 / log \x - y\dvl{y) - I log \x - y\dvl{y) - -^V2{x 



= 0, for X G supp(f2), 
< 0, for X G (-00,0]. 



(2.26) 



Since the usual equilibrium problems provide a powerful tool in the asymptotic study 
of orthogonal polynomials (cf. [8l[9l[T9]), we hope the vector equilibrium problem stated 
above will be helpful in further investigation of the asymptotics of Pk in (|2.5p : see also 

113] for a recent applications of vector equilibrium problems in some random models. 

Finally, we give the explicit formulas of the external fields Vi and V2 defined in ()2.24p , 
and the densities of the measures i^^, i^l given in ()2.16p and ()2.18p . respectively. 

Theorem 2.5. For every p,q > 0, we have 



Vi{x) = \/ q"^ + ix — plog(4x) — q\og{\/ q^ + ix + q) 

-2p- g +plog(4p2 + Apq) + glog(2p + 2q), (2.27) 
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Figure 2: The densities of the measures and z^^- (C=l) P=l-7, ^=8) 



and 



V2ix) 



0, for X < —q'^/4:, 

-\/g^+4a: 



T/ie densities of the measures u^, v\ are given by 

dv\ , . (zi + fx, s) — 2i_(x, s)) „ , , . 

X = c^^^V^V^T^V^' /orxesupp , 2.29 
ax zi+^x, sj2i_(x, sj 

dj4 ^ fc^Sfe^Sgf, e ^upp(4) and - < X < 0, 

' b^tifefeSf + 4P> forx G supp(4) and x < -gV4, 

(2.30) 



where 



^_ i{i + + p + q) 



Figure [2] shows the graph of the densities of v\ and vj^ 



Remark 1. If we take p,q ^ 0, which corresponds to fixed parameters a and the 
symbol (|2.14|) becomes 

„2 

^ - 

z 

An straightforward calculation using (j2.29p gives 



z =z + 3s^ + + 



^^a^ I 0, elsewhere, 
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Figure 3: The densities of the measures and 1^2 in the case p = q = 0. 
with 

, , , 3V3 (1 + yi^Y/^ - (1 - yi^y/' 

^^y^ = ' 

which agrees with Theorem 2.7 in [4J. This case is ihustrated in Figure [3l 

The rest of this paper is organized as follows. We first prove Proposition 12.31 in 
Section [3l The proof of Theorem 12.41 is given in Section HI We conclude this paper 
with the proof of Theorem \2.5\ where we use a nonlinear transformation to evaluate the 
integrals used to define the external fields and the equilibrium vector. 



3 Proof of Proposition 12.3 



The symbol (|2.14|) with the functions 6(s), c(s) and d{s) from (|2.13|) allows for a 
factorization 

^ ^ {z + s{s+p)){z + s{s+p + q)){z + {s + p){s + p + q)) 

By p.ip . it follows that Ag has three negative simple zeros ri,r2,r3. We order them so 
that 

ri < r2 < rs < 0; 

see Figure m for the graph of As{z). 
The derivative of As{z), 

A'^(z) = 1 - c(s)z~2 - 2(^(s)z"^ 

has three roots in the complex plane. From Figure [H we see that all zeros of A'^ are real. 
We denote the zeros of A'{z) by yi, 1/2 and 7/3 so that 

yi < y2 < < ys, (3.2) 
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as indicated in Figured! To emphasize the dependence on s, we also write yi(s), y2is) 
and y3(s). 



AM 


j 


















ri(s) 






- H 






/ i \ ^ / 








Ai yi r2\^^/r3 




y-i 


► 




- -■ 















Figure 4: The graph of As{z) ([LB]) , z G M. 



Before proving Proposition 12 . 31 we first need two lemmas. The fact that all the zeros 
Vj of the symbol Ag are strictly negative plays a key role in these proofs. 

Lemma 3.1. Assume that zi,Z2 G C are such that z\ ^ Z2, \zi\ = \z2\ and As{zi) = 
As{z2) = X. Then zi = Z2 and x G M. 

Proof. This lemma is essentially Lemma 4.1 in [14] . For the convenience of the readers, 
we repeat the proof here. 

Since zi ^ Z2 and | zi | = | Z2 1 , we may assume zi = pe^^^ , Z2 = pe^^^ with p > and 
Si / O2, 6*1^2 G ["'''"i'^]- Since the zeros rj of As are strictly negative, it is easily seen 
that f{6) := \As{pe'^^)\ is an even function on [— 7r,7r], which is strictly decreasing as 9 
increases from to vr. Hence, the equality 

\As{pe^'^)\ = \As{pe^'^^)\ 

holds if and only if 02 = —9i- It then follows that ^1 = 22, and 

X = As{zi) = As{zi) = As{z2) = X, 

so that X €R. □ 

Lemma 3.2. For each s > 0, we have ri(s) U r2(s) C M and Ti{s) fl r2(s) = 0. 

Proof. To show ri(s)ur2(s) C M, we consider two cases, based on whether the equation 
As{z)—x = has a double root or not. If x G ri(s)ur2(s) and As{z) — x = has a double 
root, then there exists zi G C such that As{zi) = x and A'g{zi) = 0. Since all zeros of 
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A'g(z) are real, it follows that x G M. On the other hand, suppose x £ Ti{s) U T2{s) and 
As{z) — X = does not have a double root, then there exist zi, Z2 £ C such that zi ^ Z2, 
\zi\ = \z2\ and x = As{zi) = As{z2). We then conclude form Lemma l3.ll that x E M. 
This proves that Ti{s) U r2(s) C M. 

To show Ti{s) n T2{s) = 0, we observe that, if x G ^lis) n r2(s), there exist three 
solutions of As{z) = x, one negative solution zi < and two complex conjugated 
solutions Z2 and Z2 ■ Moreover zi ^ Z2 and zi ^ Z2. Now \zi\ = \z2\ {x £ ^i{s)) and 
As{zi) = As(z2)- Again, by Lemma l3.H we obtain zi = Z2, which is a contradiction. 
Therefore, Ti{s) n T2{s) = 0. □ 



Proof of Proposition 12.31 Since p and q are positive, there are three local extrema 
of As{z), namely 7(5), 77(5), such that 

r?(s) < < (3{s) < 7(5). 

If X G (r/(s), /3(s)) U (7(5), 00), there exist three different real solutions of As{z) = x. 
These solutions differ in absolute value, which is obvious if x E (77(5), /3(s)) (cf. Figure 
H]) and a consequence of Lemma [3. II if x G (7('S), 00). On the other hand, there is one real 
and two complex conjugated solutions whenever x G (—00, r/(s)] U [/3(s), 7(5)]. Therefore 

Tiis) U r2(s) C (-00, rj{s)] U [/3(s), 7(s)]. 

Note that As{z) = 7(s) has a double root at 2/3 = 1/3(5) > and one negative root whose 
absolute value is less than 1/3(5). Thus 7(5) G ri(5). By the same argument, we see 
(3{s) G ri(5). In view of the fact that ri(5) is connected (see [ll],[3l Theorem 11.19]), 
it then follows that ri(5) = [/3(5), 7(5)]. Similarly, we notice that As{z) = r]{s) has a 
double root at 2/2(5) < and a negative root whose absolute value is larger than |2/2(5)|. 
Therefore r]2{s) G r2(s) and r2(5) = (—00,7/(5)]. 

To show that 7(5) is an increasing function, we introduce 

B{z, 5) = Aj'^'+P^('+P + 'ih = z{z+p){z + p + q)^ ^33^ 



(3.4) 



5 

Taking the partial derivative of B{z, s) with respect to 5, we obtain 

dB(z,s) z{z + p){z + p + q) 
ds {z — s)2 



As a function of z, it is easily seen that B{z, s) has a local minimum at 5 + 5(5 + p){s + 
P + Q)/y3{s) and 

j{s) = Asivsis)) = B{s + s{s + p){s + p + q)/y3{s),s). 
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This, together with (|3.4|) . imphes that 

dB{s + s{s + p){s + p + q)/y3{s), s) d{s + s{s + p){s + p + q)/y3is)) 



dz ds 
_^ dB{s + s{s + p){s + p + q)/y3is), s) 
ds 

dB{s + s{s + p){s + p + q)/y3{s), s) 
ds 

Z3{s){h{s) +P){z3{s) + p + q) 
{h{s)-sY 



(3.5) 



where -23(5) = s + s{s + p){s +p + q) /y3{s). As y3(s) > and p,q > 0, we have -23(5) > 
for ah s > 0. Therefore 7(5) is increasing on (0, cxd) by ()3.5p . The monotonicity of f3{s) 
and ri{s) can be proved in similar manners. Indeed, by the same argument, we have 

^ = (his) - sr ' ^'-'^ 

n,( ^ _ ziis)iziis) +p)izi{s) +P + q) , . 

{zi{s)-sr ' ^ ^ 

where Zj{s) = s + s{s + p){s + p + q)/yj{s), j = 1,2. Note that y2{s) and yi(s) are focal 
extreme points of As{z), whose zeros are —s{s+p), —s{s+p + q) and —{s+p){s+p + q). 
Therefore, in view of ()3.2p (see also Figure H]), we have 

-s{s + p + q) < y2{s) < -s{s + p) and - {s + p){s + p + q) < yi{s) < -s{s + p + q), 
which implies 

— {p + q) < Z2{s) < —p and — p < zi{s) < 0. (3.8) 

Combining ()3.6p - ()3.8p . it follows that r]'{s) > and /3'(s) < for all s > 0, which gives 
the desired monotonicity of /3(s) and r]{s). 

Finally, we come to the boundary values of 7(s), /3(s) and ??(s). A straightforward 
calculation yields that yi{s), y2{s) and ^3(5) have the following behavior as s — )• 



yi(s) = -^p{p + q){2p + q)s + 0{s), 

2p{p + q) 2x 
2/2(s) = — — s + 0{s ), 
2p + q 

ysis) = yjp{p + q){2p + q)s + 0{s). 



Hence, 



7(s) = AMs)) = p{p + q) + 0{s), 
,^{s)=AMs)) = -\q^ + 0{s), 
(i{s)=AMs))=p{p + q) + 0{s), 
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as s — 7- 0+. This proves the hmits 



hm p(s) = hm 7(5) = p(p + q) and hm r/(s) 
s— >0+ s->0+ s— 5>0+ 



On the other hand, note that yi{s), 2/2 (-s) and 2/3(5) have the following behavior as 
s —7- 00, 

yi{s) = -s^ + cis + 0{l), 
y2{s) = -s^ + C2S + 0{l), 

y,{s) = 2s' + ^P^s + 0{l), 



where ci = (— 4p — 2q — y^p"^ + pq + ^^)/3 and C2 = (— 4p — 2(jf + ■\/p'^ + pq + 9^)/3, the 
limits of 7(5), /3(s) and r/(s) as s — )• 00 can be obtained in a manner similar to the 
situation s — t- 0+, we omit the details and this completes the proof of the Proposition 

□ 



4 Proof of Theorem 12.4 



From the definitions of uf and i^l (|2.15p and (|2.18p . it is clear that supp(z^^) C 



A — 1 ^r^A gupp(jy| 

that the sets ri(s) = supp(^f) and r2(s 
Proposition 12. 3p . it follows that 



[0,oo), / dul = 1 and supp(z^2) ^ ("00,0], / di'2 = 1/2. Indeed, on account of the fact 

supp(/i2) are increasing as s increases (see 

(4.1) 



supp(z/f) = Uri(s) = ri(0 



and 



supp(i/|) 



(4.2) 



Thus, in order to show that (z^^, ^l) is the minimizer of the energy functional (|2.23p 
under the conditions stated in Theorem 12.41 it suffices to prove that the vector {vf,i>2) 
satisfies the variational conditions (|2.25p and (j2.26p . 

The basic idea is, as mentioned in Subsection 12.21 to integrate the variational con- 
ditions ()2.20p and ()2.2ip with respect to s from to ^. Here, we need a more general 
expression for the variational conditions (|2.2Up and (|2.2ip . namely 

zi{x,s) 



2 / log \x - yW\{y) - J log \x - yWl{y) -f = log 
2 / log \x - y\d^il{y) - / log \x - y\dn\{y) = log 



Z2{x,s) 

Z2{X,S) 



zz{x,s) 



(4.3) 
(4.4) 



for all 2; G C, which are contained in the proof of Theorem 2.3 of [lOj. These conditions 
reduce to (|2.20p and ()2.2ip whenever x G ^i{s) and x G r2(s), respectively. 
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Proofs of (|2.25|) and (|2.26|) . Multiplying both sides of (|4.3p by 1/^ and integrating 
with respect to s from to ^, we obtain from interchanging the order of integration that 



2 j log \x - y\duf{y) - j log \x - y\du^{y) -£=^ j log 



zi{x,s) 



Z2{X,S) 



ds. 



(4.5) 



for all X G C and some constant £ G M. The measures and t'l (|4.5p are defined in 
(|2J6]1 and ([238]) . respectively. 

Let X > 0. Since |2:i(x,s)| > |z2(a:^,s)| for every s, it follows from ()4.5p and ()2.24p 
that 

zi{x,s) 



2 J \og\x - y\dvl{y) - j \og\x - y\dvl{y) - 1 <^ j log 



Z2ix,s) 



ds = \vi{x). (4.6) 



Suppose now x E supp(i/^), then x G ri(^) by (|4.ip . and therefore x G ^i{s) for every 
s > ^, since the sets are increasing. Thus |zi(x,s)| = |2;2(x,s)| for every s > ^, and 
equality holds in (|4.6p for x G supp(i/^). This completes the proof of ()2.25p . 

The variational condition ()2.26p for z^l can be proved in a manner similar to (|2.25p 
by using (j4.4p and (|4.2p . we omit the details. □ 



5 Proof of Theorem 12.51 

We conclude this paper with the proof of Theorem 12.51 We start with the following 
lemma that embodies the behavior of three solutions of As{z) as s — t- 0+. 

Lemma 5.1. Let As{z) he given by (|3.ip . and let zi{x,s), Z2{x,s) and Z3{x,s) be the 
solutions of As{z) = x, ordered as in (jl.Sp . Then 



lim zi(x, s) = X — pip + q), 
s-i>0+ 



lim s ^Z2{x,s) 



p{p + q){2p + q + y/(f+Ax) 

2(^p2 _|_ pg, _ 3.^ 



lim s ^zs(x,s) 



p{p + q){2p + q - \Jq^ ^ 4x) 



2(p2 j^-pq — x) 

Proof. The lemma follows by a straightforward computation. 

To prove Theorem 12.51 we need to establish the identities ()2.27p - ()2.30p . 



(5.1) 



□ 



Proof of ()2.27p . Let x > 0. From Proposition 12. 3|, it follows that there exists a unique 
s*(x) > so that for all s > 0, 



X £ ri(s) <^=^ s > s*{x). 
Then \og\zi{x, s)/ Z2{x, s)\ = for ah s > s*{x), and so by ()2.24p 



Vi{x) 



s*{x) 



log 



zi{x, s) 

Z2{X, S) 



ds. 



(5.2) 



(5.3) 
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There is a special value 

^0 = pip + q) = li™ /3(s) = lim 7(5) 

that belongs to every ri(s) for any s > 0. Then s*(xo) = and 

Vi{xo) = 0. (5.4) 

The derivative of (15.311 is 



VT(., . /'■'" (-1-^ - -i^^) (5.5) 

7o \zi{x,s) dx Z2{x,s) dx J 
In order to evaluate this integral, we introduce new variables 

Zj{x,s) = + s, J = 1,2,3. (5.6) 

Zj[x,s) 

Since each Zj{x,s) is a solution of As{z) = x, it follows that Zj{x,s) for j = 1,2,3 is a 
solution of the equation B{z, s) = x, where 

Biz,s) = '-^^±^^^^±P±^; (5.7) 
z — s 

see ([3:3]) . 

Taking partial derivatives with respect to s and x on both sides of B{zj{x, s), s) = x, 
and applying the chain rule, we obtain 



dB , . ,\ dzj{x,s) dB . , , 
— {zj{x,s),s) j + — (zj (x,s),s) = 0, 



for j = 1, 2, 3. From (|5.7|) . it is elementary to deduce that 

— iz ix s) s) = + + P + ^ ^ 

ds '' ' {zj{x,s) — sy Zj{x,s) — s' 

Combining this with (j5.8p . we obtain 

X dzj (x, s) (9ij' (x, s) 



(5.8) 



Zj{x, s) — s dx ds 
By ()5.6p . we also note that 



j = 1,2,3. (5.9) 



dzj{x,s) s{s + p){s + p + q) dzj{x, s) Zj{x, s) — s dzj{x, s) 
dx Zj{x,s) dx Zj{x,s) ds 



(5.10) 
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Substituting (|5.1U|) into (|5.9|) gives 

1 dzj{x,s) ldzj{x,s) 



Zj{x,s) ds X ds 

Hence, using (jS.lip in (|5.5p we get 



j = 1,2,3. (5.11) 



ds ds 

It then follows from the fundamental theorem of calculus that 

V((x) = -[zi(x,s*(x)) -Z2(x,s*(x)) - lim (zi(x,s) -Z2(x,s))]. (5-12) 

X s^0+ 

By definition, s*{x) is the smallest value of s > for which x G ri(s). Then 
X = 7(s*(x)) if xo < X and x = /3(s*(x)) if < x < xq. We can observe from Figure S] 
that 21(7(5), s) = 22(7(5)55) and zi(/3(s),s) = 22(/3(s),s). Therefore, 

i'i(x,s*(x)) = Z2ix,s*{x)) 

and (|5.12p reduces to 

VUx) = lim (zi(x, s) — 'z2(x, s)). (5.13) 

X s-^0+ 

From (|5.6p and Lemma |5.H we find that 

2(p^ + pq — x) 



lim zi(x,s)=0 and lim Z2{x,s) 



s^o+ s^o+ 2p + q+ y/^~+4x 

Then (|5.13p leads to 

V(ix) = -1 '(P'+P^--^ . (5.14) 
X2p + q + V g2 + 4x 

We obtain Vi{x) by integrating (j5.14p with respect to x, which gives 



Vi{x) = y/q'^ + 4x -plog(4x) - qlog{y/ q"^ + 4x + q) + C. 

The constant of integration C can be determined by requiring Vi{xo) = 0; see ([5 
This gives 

C = -2p - q + p log(4p^ + Apq) + q log(2p + 2q), 
and ([2:271) is proved. □ 
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Proof of (|2.28|) . Let x < 0. Again, it follows from Proposition 12.31 that there exists a 
unique s*{x) > so that for all s > 0, 



X E r2(s) <^=^ s > s*{x). 
Then log 1 2:2(3;, s)/2;3(x, s)| = for ah s > s*{x), and so by (|2.24|) 



V2ix) -- 

There also exists a special value 

Xo 



s*{x) 



log 



Z2{X,S) 



Z3{X,S) 



ds. 



(5.15) 



(5.16) 



-q /4 = lim rjis) 



so that X G ^2{s) for any s > if x < xq. Then, for any x < xq, we have s*{x) = and 

V2{x) = 0. (5.17) 

If — = Xq < X < 0, by the same argument in the proof of (j2.27p . we have 

1 



1^2(2;) = -- lim (2;2(x, s) - zs{x, s)), (5.18) 

X s^0+ 

where Zj{x, s), j = 2, 3 are given in (|5.6p . From Lemma |5. 11 we see 

2(jP'+pq — x) _ 2{p^+pq — x) 



lim Z2{x,s) = - 

«^o+ 2p + q + yjq^ + 4x 



and lim 2:3 (x, s) = - 

2p + q- V g2 _^ 43, 



Hence, ()5.18p leads to 



\Jq^ + 4x 



(5.19) 



By integrating (j5.19p with respect to x, we obtain 



V2{x) = -2./JT^x + q\og + 

^9 - VT + 4x^ 

for — (?^/4 < X < 0. The constant of integration C can be determined by requiring 
V2{xq) = 0; see (fCTTl) . This leads to C = and (f?:^ is proved. □ 

Proofs of ([2:29]) and (|230]) . By the definitions of i/f and i/| in (|2T6]) and (|238]1 . we 
obtain from (fTTO]) and (f2T7P that 



dx ^ ' 2TTi^ 



1 dzi-{x,s) 



1 dziJ^{x,s) 



1 /■« 



zi_(x,s) 9x 2;i+(x,s) dx 

1 5z2+(x,s) 1 9z2-(2;,s) 



2vri^ Jo \-2^2+(a;;s) 9x 2;2_(x,s) 5x 



ds. 



ds. 
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With s*{x) given in (|5.2p and (|5.15|) . it is readily seen that 



^^1 _ 1 ( 1 dzi-{x,s) _ 1 dzi+{x,s) , 



dx liri^ J gtf^^-j \zi-{x, s) dx zi+{x,s) dx 

dyj / N _ 1 /"^ / 1 g2:2+(x,5) 1 gz2-(a^,g) \ 

dx 2TTi^ J \z2+{x, s) dx Z2-{x,s) dx 



Now, by introducing the change of variable (|5.6p . one can evaluate the above integrals 
using similar methods as given in the proofs of (|2.27|) and (|2.28|) . We omit the details 
and this completes the proof of Theorem 12.51 □ 
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